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SOME INEQUALITIES FOR THE TRIGAMMA FUNCTION IN 
TERMS OF THE DIGAMMA FUNCTION 

FENG QI AND CRISTINEL MORTICI 


Abstract. In the paper, the authors establish three kinds of double inequal¬ 
ities for the trigamma function in terms of the exponential function to powers 
of the digamma function. These newly established inequalities extend some 
known results. The method in the paper utilizes some facts from the asymp¬ 
totic theory and is a natural way to solve problems for approximating some 
quantities for large values of the variable. 


1. Motivations and main results 
The classical Euler gamma function may be defined for 5i(2) > 0 by 

/> OO 

r(*)= / t*- 1 e"*dt. 

Jo 

The logarithmic derivative of the gamma function r(^) is denoted by 

V’O) = -r-[inrO)] = yy-y 
dz r(z) 

and called the digamma function. The derivatives ip'(z) and ip"{z) are called the 
trigamma and tetragamma functions respectively. As a whole, the functions (z) 
for k G {0}UN are called the polygamma functions. These functions are widely used 
in theoretical and practical problems in all branches of mathematical science. Con¬ 
sequently, many mathematicians were preoccupied to establish new results about 
the gamma function, polygamma functions, and other related functions. 


1.1. The first main result. In 2007, Alzer and Batir [3, Corollary] discovered 
that the double inequality 




x exp 


—x — —^(x + a) 


< r(®) < V^Trx x exp 


-x- -ip{x + P) 


( 1 . 1 ) 


holds for x > 0 if and only if a > I and fi < 0. For information on generalizations 
of results in the paper [3], please refer to [12, 32] and references cited therein. 
Motivated by the double inequality (1.1), Mortici [27] proposed the asymptotic 
formula 


r(x) ~ vTre b (x + b) x exp 



1 


ip(x + c) 


x —y oo 


and then determined the optimal values of parameters 6, c in such a way that this 
convergence is the fastest possible. 
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In 2004, Batir [4, Theorem 2.1] presented the double inequality 


T(c) exp^aOe 1 ^ - e^ x > + 1 ] < T(*) 

< T(c) exp| + l]|, ( 1 . 2 ) 

for every x > c, where c = 1.461... is the unique positive zero of the digamma 
function ip and 7 = 0.577... is the Euler-Mascheroni constant. In 2010, Mortici [25] 
discovered the asymptotic formula 


r(x) ~ v/ ^' exp [ip(x)e^^ — e+ l], x —> 00 . 

In 2011, the double inequality (1.2) was generalized by [11, Theorem 2] to a mono¬ 
tonicity property which reads that the function 


fs,t( X ) 


_ g s ,t ( x ) _ 

WsA x ) ~~ !] exp[3( jt (a')] + 1’ 

1 


x ^ c 


g'Ucy 


X = c 


on ig (—a, 00 ) is decreasing for |f — s| < 1 and increasing for \t — s| > 1 , where s 
and t are real numbers, a = min{s, t}, c € (—a, 00 ), and 


9s,t{ x ) 



In 


T(m + t) T(c + s) 
T(u + s) T(c + t) 


u + s) — il>(c + s)] d u, 


dtt, 


s y t 
s = t 


onig (—a, 00 ). 

In 2000, Elezovic, Giordano, and Pecaric [ 8 ] found the single-sided inequality 


ip'(x ) < e ^ x \ x > 0 . 


(1.3) 


This inequality is closely related to the monotonicity and convexity of the function 

Q(x) = e ^ x+1 "> - x 

on (—l,oo). See also [13] and plenty of references therein. By the way, as a 
conjecture posed in [15, Remark 3.6], the complete monotonicity of the function 
<2( x) on (0, 00 ) still keeps open. An infinitely differentiable function / is said 
to be completely monotonic on an interval I if it satisfies (—1 ) k f^ k \x) > 0 on 
I for all k > 0. This class of functions has applications in the approximation 
theory, asymptotic analysis, probability, integral transforms, and the like. See [7, 
Chapter 14], [22, Chapter XIII], [51, Chapter 1], and [52, Chapter IV], For more 
information on the functions 

Q'(x — 1) = 'll)'{x)e^ ( ' x ^ > — 1 and Q"(x — 1) = {ip"(x) + [^(a;)] 2 }^^, 

please see the papers [14, 15, 17, 33, 36, 37, 38, 41, 42, 43, 48, 44, 47, 54, 55], the 
expository and survey articles [34, 35, 45, 46], and a number of references cited 
therein. In 2011, Batir [5, Theorem 2.7] obtained the double inequality 

(x + a*)e~ 2 ^ x+1) < iP'{x + !)<(* + b*)e~ 2 ^ x+1 \ 


x > 0 , 


(1.4) 
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where the constants a* = j and b* = = 0.518... are the best possible. In 

other words, Batir [5] proposed the approximation formula 

i/(x + l) ~ (x + a)e~ 2 ^ x+1) , (1.5) 

where a is a constant. A numerical computation shows that the approximation (1.5) 
gives a better result when choosing a = a*, rather than the value a = b*. This 
fact is somewhat expected as Batir obtained (1.4) as a result of the decreasing 
monotonicity of the function 

9{x) = ip'(x + l)e 2 ^^ +1 ^ — x 

on (0, oo), with 0(0) = ^77 and 9{ 00 ) = i. Hence, the function 9(x) becomes 
gradually closer to 0(oo), as x approaches infinity. The decreasing monotonicity of 
the function 9(x ) may be rewritten and extended as that the function 

Q{x)='il)\x)e^ {x) -x + ^ 

is decreasing and positive, but not convex or concave, on (0,oo). In 2010, among 
other things, Qi and Guo [43, Corollary 2] proved that the function 

f P , q (*) = e^* +1 ) - qx 

for p ^ 0 and q £ R is strictly convex with respect to r € (— 1 , 00 ) if and only if 
p > 1 or p < 0. Since 

] p aip(x) j2 aip(x) 

-^- = atp’(x)e a ^^ and --— 5 —= a{ip"(x) + a[il)'(x)] 2 }e a ^^ x \ ( 1 . 6 ) 

d x cl ar 

the above functions between (1.3) and (1.6) have something to do with the more 
general function 

fa, b,a,f),\(x) = e ^ iX+b) +aX 2 +pX + \ 

and its monotonicity and convexity, even its complete monotonicity. 

For obtaining accurate approximations of the type (1.5), the value <1 = \ should 
be used. However, approximations of the form 

i/j'(x + 1 ) ~ [x + a{x)]e~ 2 ^ x+1 \ 

with a(x) —> \ as x tends to 00 , are much better. Furthermore, several experiments 
using some asymptotic expansions lead us to the claim that approximations of the 
type 

^'{x + 1 ) ~ [x + a(z)] exp -2ip(x + 1 ) - 
are more accurate. Consequently, we obtain the following theorem. 

Theorem 1.1. For x > 3, the double inequality 


[x + a{x)\ exp —2ij)(x + 1) — 


< if'(x + 1) 


< [x + /3(x)\ exp —2ij)(x + 1) — 


is valid, where 


. , 1 1 
“W = o + bb- 


90a; 3 60a; 4 


and j3{x) = — + 
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1.2. The second main result. In 2013, Guo and Qi [14, Lennna 2] proved that 

ip\x) < e 1/x - 1 , x > 0 . ( 1 . 8 ) 

This inequality has been generalized to the complete nronotonicity of the function 

g 1 /* _ ^'(x) — 1, x > 0 (1.9) 

and others. For detailed information, please refer to [16, 21, 39], [40, Theorem 3.1], 
[49, Theorem 1.1], [50, Theorem 1.1], and closely related references therein. In some 
of these references, it was obtained that 

lim [e 1 /® — ip'(x )] = 1 . 

x—too L J 

The approximation ip'i x ) ~ e 1 ^ ~ 1 indeed gives a good result for the large value 
of x. Now we propose the improvement 

i/j'{x) ~ e 1/x+ ^ (x ) - 1 , 

where p(x) —>• 0 and p(x) = o(x) as x —> oo. Our main result may be stated in 
details as the following theorem. 


Theorem 1.2. For x > 3, we have 

e m(x) _ x < < e M(x) _ 1; (!. 10 ) 


where 


m(x) 


1 1 

x 24a; 4 


7 

360a; 6 


and M (x) 


m(x) + 


1 

90a; 7 ' 


1.3. The third main result. In 2014, Yang, Cliu, and Tao found [53, Theorem 1] 
that the constants p = 1 and q = 2 are the best possible real parameters such that 
the double inequality 

0(x,p) < ip'{x + 1 ) < 0(x,q) ( 1 . 11 ) 

holds for x > 0 , where 


9(x, m) 


e m/(x+l) _ g -m/x 

2 in 


Because 

ip\x+l) = ip'(x) - -1, 

X A 

the double inequality ( 1 . 11 ) may be reformulated as 


( 1 . 12 ) 


+ 0(x, 1) < ip'(x) < —j + 9(x, 2), x > 0. (1.13) 

ar x z 

The complete nronotonicity of the function (1.9) implies that 

e — ip'{ 1 ) > e 1 ^ x+1 - ) — ip'{x + 1 ) > 1 , x > 0 


which may be rearranged as 

gi/O+i) _ e ^ ip'(l) < tp'(x + 1 ) < e 1 ^ x+1) — 1 < - sinh — (1.14) 

2 x 

on (0, oo). When t < 1.6 ..., the lower bound e l ^ x+1 ' > —e + 'ip'il) is better than the 
corresponding one in [53, Corollary 2]. The upper bound — 1 in (1.14) and 

the upper bound 9(x,2) in (1.11) can not be compared with each other. However, 
the upper bound + 9(x,2) in (1.13) is better than the upper bound e 1 ^ — 1 
in ( 1 . 8 ). 

In this paper, we will improve the double inequality (1.11) as follows. 
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Theorem 1.3. For x > 1, we have 


e[x ' 11 + 24? - S? < ' P ' (X +l)< S(X ’ 11 + 2ib 


°^ 2) -ik < *' lx + 1)<e{x - 2) -ik + 9b- 

2. Proofs of Theorems 1.1 to 1.3 

Now we start out to prove our three main results. 

In 1997, Alzer [2, Theorem 8 ] proved that for m,n > 1 the functions 


/ 1 \ i 2m 

F m {x) = In T(x + 1) — ( x + - J In x + x — - In 27r — 


2i(2i — l)x 2i_1 


/ 1 \ \ 2 ™ 1 

G n (x) = — lnr(a; + 1 ) + ( x + — j In a; — x + — ln27r + 


2i(2i — l)x 2t 1 


are completely monotonic on (0, oo), where Bj are the Bernoulli numbers given by 
the generating function 

t ™ & 

e*-l 2 + ^ 2j (2j)V 

j=o J j=o 1 

See also [19, Theorem 2], [20, Theorem 2.1], and [31, Theorem 3.1]. From F[(x) < 0, 
G' 2 (x) < 0, F"(x) > 0, F 2 (x) > 0, G 2 {x) > 0, and G'f(x) > 0, it follows that 

lXlX+ 2x~ 12 X 2 + 240x I ~ 252x6 <^ x + 1 ) <lnx+ 2 ^~ u x 2 + 240a; 4 ’ 

^'(x + l) > + (2.2) 

1 1 _L_ 1_ _L_ 1 _ , 

x 2x 2 6 a ; 3 30a; 5 42a; 7 30a; 9 < X 

1 - 1 1 1 1 

< x 2x 2 6 a ; 3 30a; 5 42a; 7 ’ 

and, also by ( 1 . 12 ), 

1 1 _L_ l_ J_ _1_ , 

x 2a ; 2 6 a ; 3 30a; 5 42a; 7 30a; 9 < X 

111 1 1 1 5 

< x 2a ; 2 6 a ; 3 30a; 5 42a; 7 30a; 9 66 a; 11 ' 

Proof of Theorem 1.1. By taking the logarithm, the double inequality (1.7) may be 
rearranged as 

F{ x) = In [a; + a(x)\ — 2 ip(x + 1) — Intp'fx + 1)-— r < 0 

120 a ; 4 


G{x) = - In [a; + fi(x)\ + ln^a; + 1) + 2if>{x + 1) + 4 < 0 
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for x > 3. By virtue of inequalities (2.1) and (2.3), one may deduce that F(x) < 
Fi(x) and G(x) < Gi(x), where 


Fi (x) = In [a; + a{x)\ — 2 I In x + - 


and 


1 


1 


2x 12x 2 240a: 4 252a; 6 


1 


1 


1 


1 


__ 

x 2x 2 + 6 a ; 3 30a: 5 42a; 7 30a; 9 


1 


G\ (x) = — ln[x + /3(x)\ + In (- 


1 


1 


+ 


120 a ; 4 

1 


2x 2 6 a : 3 30a: 5 42a; 7 


+ 21 In x + J- - 


1 


Since 


and 


where 


f;(x) = 


2x 12a ; 2 240a: 4 

A( x — 3) 


1 


120 a; 4 ' 


G'i{x) = 


105a: 7 (180x 5 + 90a; 4 + 2a: - 3)B(x - 3) 
C(x — 3) 


15a; 5 (90a; 4 + 45a ; 3 + 1 )D(x - 3) ’ 


A{ x) = 84000x 13 + 3753750a; 12 + 75109720.T 11 + 896904120a; 10 

+ 7160223140a; 9 + 40457327085a; 8 + 166700796732a; 7 + 507517074474a; 6 
+ 1141703970759a; 5 + 1873185114060a; 4 + 2175691772642a; 3 
+ 1690072075536a; 2 + 783944661553a; + 162974708124, 

B{ x) = 21 Ox 8 + 4935x 7 + 50750a; 6 + 298305a; 5 + 1096193a; 4 

+ 2578821a; 3 + 3792857a; 2 + 3188649a; + 1173161, 

C(x) = 23625a; 9 + 604200a; 8 + 6818475a: 7 + 44510545a; 6 
+ 184933335a: 5 + 506070905a; 4 + 909421185a; 3 
+ 1030441127a; 2 + 663679092a; + 183168418, 

D{x) = 210a ; 6 + 3675a; 5 + 26810a: 4 + 104370a: 3 + 228683a; 2 + 267393a: + 130352 


are polynomials with positive coefficients, the functions F\(x) and G\(x) are strictly 
increasing on [3, oo). Furthermore, since 


lim F\(x) = lim G\{x) = 0, 

x—too x—yoo 


it follows that Fi(x) < 0 and G\(x) < 0 for x > 3. As a result, we have F(x) 
Fi(x) < 0 and G(x) < G\(x) < 0. The proof of Theorem 1.2 is complete. 


Proof of Theorem 1.2. The inequality (1.10) can be written as 

m(x) < ln[l +ip'(x)\ < M(x). 

Considering the double inecjuality (2.4), since m(x) — ln[l + ij>'(x)\ < m\(x) and 
M(x) — ln[l + ip'{x)\ > Mi (a;), one may see that it suffices to show that 

1 11 1 

2x 2 " 


mi(x) = m(x) — ln( 1 H- 1 - 


+ 


6 x 3 30x 5 42x 7 30x 9 


< 0 


□ A 
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and 

MlW = - ‘"O + \ + W + 675 - 35 ? 

A straightforward computation gives 

° araJ MiM 

where 


1 

42a 7 


1 5 \ 

30a 9 66 a: 11 / 


> 0 . 


F(x- 3) 
l80x 8 H(x - 3) ’ 


E(x) = 980a: 8 + 22485a; 7 + 221130a; 6 + 1212855a: 5 + 4032099a: 4 

+ 8229303a: 3 + 9865371a; 2 + 6074127a; + 1290163, 

F(x) = 66495a: 10 + 2146155a: 9 + 31007240a: 8 + 263913573a; 7 
+ 1464790565a: 6 + 5537745108a: 5 + 14437981040a: 4 
+ 25626153678a; 3 + 29624987873a; 2 + 20135221233a; + 6106987838, 
G(x) = 210 a ; 9 + 5880a 8 + 73185a: 7 + 531440a: 6 + 2481255a: 5 + 7724423a; 4 

+ 16033731a; 8 + 21398207a 2 + 16660569a + 5765861, 

H{x) = 2310a 11 + 78540a 10 + 1213905a 9 + 11258170a 8 + 69614160a 7 
+ 301344043a 6 + 931827204a 5 + 2058324400a 4 + 3182887290a 3 
+ 3281444518a 2 + 2029943157a + 570820414 


are polynomials with positive coefficients. This means that rn[ (a) > 0 and M[ (a) < 
0 for a > 3, that is, the function mi (a) is strictly increasing and the function M\(x) 
is strictly decreasing on [3, oo), with the limits 

lim mi (a) = lim Mi (a) = 0 . 

tc->-oo x —>oo 

Accordingly, one obtain that mi (a) < 0 and M\{x) > 0 on [3,oo). The proof of 
Theorem 1.2 is complete. □ 


Proof of Theorem 1.3. The left hand side inequality in (1.15) can be written as 
c(a) 4 I e 1 /( x + 1 ) _ I e -V- + 1 5 

Since 


24a 5 48a 6 


— tp'(x + 1 ) < 0 . 


k—0 


7 if i x fc 


k\ V x 


and the inequality ( 2 . 2 ) is valid, we have 




where 


P(x) = 5040a 7 + 5040a 6 - 2520a 5 + 840a 4 + 210a 3 - 798a 2 + 1057a - 1. 
Hence, it is sufficient to prove that 

ci(a) = ——r - In , P ^\ < 0 . 


a +1 


5040a 7 
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A straightforward differentiation yields 

, 7630a; 2 + 6329a; - 7 _ 7630(a - l) 2 + 21589(a - 1) + 13952 

Cl X a(a + l) 2 P(a) x(x + l) 2 Q(x — 1) 

where 

Q{x) = 5040a: 7 + 40320a; 6 + 133560a: 5 + 240240a; 4 

+ 255570a- 3 + 161112a 2 + 56371a + 8868. 

This implies that the function Ci(a) is strictly increasing on [l,oo). Furthermore, 
because of lim^oo ci(a) = 0, it follows that ci(a) < 0 on [1, oo). The proof of the 
left hand side inequality in (1.15) is complete. 

Similarly, we may verify other inequalities in (1.15) and (1.16). For the sake 
of saving the space and shortening the length of this paper, we do not repeat the 
processes. The proof of Theorem 1.3 is complete. □ 


3. Remarks 

Finally we give several remarks on our three main results. 

Remark 3.1. The right hand side in (1.7) is better than the right hand side in (1.4), 
but the left hand side in (1.7) is not better than the left hand side in (1.4). This is 
because the inequalities 

1 


(a; + a*)e 2 V'( a: + 1 ) > [# + a{x)\ exp 

and 

[x + /3(a)] exp —2%l>(x + 1) — 
which are equivalent to 


—2ip(x + 1 ) — 


1 


120 a ; 4 


120 a ; 4 
< (x + b*)e~ 2 ^ x+1) 


and 


«M = ln(* + i + A5-Aj 


, v , , 1 1 

v{x) = ml x + - + 


- In [x + - - 


1 


120 a - 4 


< 0 


90a 3 


1 


— In ( a + 


6 e 27 


< 0 , 


(3.1) 


(3.2) 


120 a 4 

are valid for a > 1. The inequalities (3.1) and (3.2) may be verified as follows. A 
direct differentiation gives 

540a 6 + 60a 5 + 180a 4 + 90a(a - 4 ) 2 + ^a + 2 
u'{x) = - ^- > 0 


for a > ^ and 


for a > 1 , where 


y'(a) = 


30a 5 (2a + l)S(a — j^) 

_ Q(x - 1) _ 

30a 5 (90a 4 + 45a 3 + l)( 6 e 27 a + 7 r 2 ' 


> 0 


S(x) = 180a 4 + 162a 3 + ^a 2 + §a + ^ 

O OU iZo 


Q(x) = 2700(t r 2 - 3e 2j )x 8 + 21600(tt 2 - 3e 27 )a 7 + 75600(tt 2 - 3e 27 )a 6 
+ 180(840tt 2 - 2521e 27 )a 5 + 630(300 tt 2 - 901e 27 )a 4 


and 
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+ 45(3361?r 2 - 10096e 27 )a; 3 + 45(1683 tt 2 - 5044e 27 )a; 2 

+ 3(7245tt 2 - 21538e 27 )x + 2(1373 tt 2 - 4002e 27 ) 

are polynomials with positive coefficients. This means that the functions u(x) and 
v(x) are strictly increasing on [l,oo). Further, from the limits 

lim u(x) = lim v(x) = 0, 


it follows that u(x) < 0 on (0, oo) and v(x) < 0 on [1, oo). 

Remark 3.2 (The asymptotic series of i/j'( x + l)e 2 ^^ x+1 ^). Whenever an approxima¬ 
tion formula f{x) ~ g(x) is considered in the sense that the ratio lim^oo = 1 
(sometimes lim x _ ) . (x) [/(a;) — g(x)\ = 0), there is a tendency to improve it by adding 
new terms, or an entire series, of the form 

OO 

/(*) ~s(*) + X3- 

fc=l 

Such a series is called an asymptotic series and it plays a central role in the theory 
of approximation. Although such a series is often divergent, by truncating at the 
m-th term, it provides approximations of any desired accuracy 1 11 other 

words, the formula 


m / , \ 

/(*) ~ 9(x) + J2 ^ + 0 ( ) > 

fc=l 7 ' 


x —> 00 


is valid for every integer m > 1. 

It is well known [1] that the digannna and trigamma functions admit respectively 
the following asymptotic expansions: 


1 °° D 1 °° O 

iP(x + l)~lnx+ — and + --5 + J2-Z^- 


(3.3) 


k—2 k -1 

The asymptotic expansions in (3.3) can be written explicitly as 

1111 


Hx + 1) ~ In a; + ^ ^ 132xl0 


+ 


and 


1 


1 


1 


Z 11 1 

X ^ x 2x 2 6a; 3 30a; 5 ' 42a; 7 30a; 9 ' 66a; 11 


+ 


In order to construct the asymptotic expansion of the function ip'(x+ l)e 2 ^ x+1 ), 
we recall from the papers [6, 10] and the monographs [9, pp. 20-21] and [18, 
pp. 539-541] the following classical results from the theory of asymptotic series: 
(1) if 

OO 

Pk 
x k 


and 


k =0 


k =0 


as x —> oo, then 


where 


u(x)v(x) ~ 


r_k_ 

rv* k ' 


(3.4) 


k—0 

Tk = ^2 Piq j’ 

i-\-j=k 
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(2) if 


then 


where ao 


x ^ oo, 


k=0 


exp/(x) ~ ^ z-too, 

/c—0 

expao and 

**=e j e n^ 

j=t J h+-+ij=ke=i 


(3.5) 


for k > 1. 

Now the asymptotic series of the function ^'(x + l)e 2 ^( a:+1 ) can be computed in 
two steps. Firstly, by virtue of the formula (3.5), we may transform the series (3.3) 
to obtain the series of e 2 ^^ +1 ^. Secondly, with the help of the formula (3.4), we 
multiply the series of e 2 ^^ +1 ^ and the second series in (3.3). 

Because the general term (in terms of the Bernoulli numbers) of the series of the 
function %j)'(x + l)e 2 ^’( a:+1 ) has an unattractive form, we just write down the first 
few terms as follows: 


il>\x + l)e 2 ^ x+1) 


1 , 1 1,2 

2 + 90x3 “ + 567x5 


43 

2268x6 



Till now we can see immediately that the functions a(x) and /3(x) in Theorem 1.1 
are truncations of the series (3.6) at the first three or four terms. 

When more terms of (3.6) are considered, the conclusion of Theorem 1.1 remains 
true and the estimates become more accurate. 


Remark 3.3. We conjecture that the function e M ^ — ^'(x) — 1 is completely mono¬ 
tonic on (0, oo), but the function ip'(x) — e m ^ + 1 is not. 

Remark 3.4. The reason why the double inequality (1.11) is of interest is because it 
provides much accurate estimates for the trigamma function ip'{x) for large values 
of x, that is, 

ip'(x + 1) ~ 0(x, to), x—>■ oo. (3.7) 

From (1.15) and (1.16), it is easily deduced that 

ip'(x + 1) - 9(x, 1 ) = 0 

tlren it is easy to see that the best approximation of the form (3.7) may be obtained 
by taking m = 2. 


^ and ip'(x + 1) - 0(x, 2) = 


Remark 3.5. Evidently, our inequalities (1.15) and (1.16) are much stronger than 
those in [53]. Our proofs for (1.15) and (1.16) are also much simpler than the 
original proof of [53, Theorem 1], This shows that the natural approach to solve 
problems of approximating some quantities for large values of the variable is the 
theory of asymptotic series. This method or approach has been utilized and applied 
in the papers [23, 24, 26, 28, 29, 30], for examples. 


Remark 3.6. In this paper, we essentially talk about the relations between inequal¬ 
ities, asymptotic approximations, and complete nronotonicity. 
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